Hidden past of dark energy cosmological models 
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In this paper we analyse the possibility of having homogeneous isotropic cosmological models with 
observers reaching t = oo in finite proper time. It is shown that just observationally-suggested dark 
energy models with w 6 ( — 5/3,-1) show this feature and that they are endowed with an exotic 
curvature singularity. Furthermore, it is shown that non-accelerated observers in these models may 
experience a duration of the universe as short as desired by increasing their linear momentum. A 
subdivision of phantom models in two families according to this behavior is suggested. 
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I. INTRODUCTION 

During the last years there has been mounting experi- 
mental evidence from different sources (supernovae type 
la [l|, redshift of distant objects [|J and temperature fluc- 
tuations of background radiation Q) supporting an ac- 
celerated expansion of our Universe at present time (cfr. 
for instance Q for a review). 

Trying to explain this fact, several proposals have been 
made, such as dark energy contents for the universe or 
modifications of the theory of gravity, which have pro- 
duced a menagerie of new types of singular events in the 
respective cosmological models, traditionally restricted 
to Big Bang and Big Crunch singularities. For instance, 
we may find in phantom energy models Big Rip singular- 
ities Q . One of this models has been shown to be stable 
against quantum corrections ffj}. An attempt to explain 
the accelerated expansion without violating all energy 
conditions [7] produces sudden singularities. Most re- 
cently, inaccessible singularities in toral cosmologies have 
been added to the list Q. 

There have been several attempts to organize these 
families of singular events in thorough classifications. In 
@ all types of singular events in Friedmann-Lemaitre- 
Robertson- Walker (FLRW) cosmological models are clas- 
sified according to the coefficients and exponents of a 
power expansion in time of the scale factor of the uni- 
verse around the event. In [l(| singularities are classified 
using the finiteness of the scale factor, the density and the 
pressure of the universe. In [Tl| the behavior of causal 
geodesies close to singular events and the strength of the 
singularities are analysed. 

This line of research has proven successful showing un- 
expected features of FLRW cosmological models near the 
singularities. For instance, it has allowed to show that 
sudden singularities are weak [T^, since tidal forces do 
not disrupt finite objects falling into them [u3. [l5j. 



Another intriguing feature concerning Big Rip singu- 
larities is that photons do not experience such fate for 
effective equations of state, p = wp, with w € (—5/3, — 1) 
(that is, those comprised between the superphantom [16[ 
and the phantom divide), since they require an infinite 
lapse of time to reach that event [ll| . Since this range of 
the parameter w comprises the observationally accepted 
values [17], which are slightly below -1, this fact is more 
than a mere curiosity. 



Following the idea of classifying the singular events 
arising in FLRW cosmological models, it is worth men- 
tioning that all classifications are incomplete in a sense: 
they unveil what happens at a finite coordinate time t, 
but they are elusive when asked about infinite t. This 
may seem a pointless consideration, since in most cases 
an infinite coordinate time lapse corresponds to an infi- 
nite time lapse experienced by the observer, but the men- 
tioned example about photons in phantom cosmologies, 
where a finite coordinate time lapse requires an infinite 
proper time shows us that the issue is far from being 
trivial. 
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To this aim in Sec. |TT] the equations governing causal 
geodesies in FLRW cosmological models are reviewed. 
In Sec. IIIII the conditions for a causal geodesic to reach 
t = ±oo in finite proper time are derived. It will be 
shown that just phantom models fulfill this property. In 
Sec. IIVI it will be discussed if this abrupt end of causal 
geodesies is an actual singularity or not. Analysis of the 
Ricci curvature as measured by the observers will settle 
the issue, in spite of the zero value of curvature scalar 
polynomials there. In fact these are strong curvature 
singularities. Finally, the consequences of these facts will 
be discussed in Sec. IV] 
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II. GEODESICS IN FLRW COSMOLOGICAL 
MODELS 

The metric for FLRW cosmological models may be 
written, 



ds 2 
f(r) 



-dt 2 +a 2 (t) {f(r)dr 2 +r 2 (d6 2 + sin 2 0d<p 2 )} 
' fc = 0,±l, (1) 



vT~— kr 2 ' 



in terms of spherical coordinates r, 9, (j) with their usual 
ranges and a coordinate time, with a range depending on 
the type of cosmological model. 

Three families of models are comprised in this expres- 
sion, open models with k = —1, flat models with k = 
and closed models with k = 1. Observations favor flat 
models, but we keep for our purposes the general formula. 

Free-falling observers in a spacetime are modeled by 
timelike geodesies parametrized by proper time r, since 
these curves have the property of vanishing acceleration. 
The use of proper time allows us to write the velocity 
u of the parametrization of the geodesic (i, r, 9, 4>) as a 
unitary vector, 



x , x J 



(2) 



where the dot means derivation with respect to t. 

There are three types of geodesies: timelike (<5 = 1), 
spacelike (<5 = —1) and lightlike (<5 = 0). We consider 
just causal geodesies, S — 0, 1, since they are the only 
ones that may carry signals or observers. 

A quick way to write down simple geodesic equa- 
tions for these spacetimes is taking into account that 
the universe is homogeneous and isotropic and therefore 
geodesies are straight lines in the spacetime and we may 
take 9 = = <p without loss of generality. 

It is also easy to check that the vector <9r = <9 r //(r), 

arcsinhr k = — 1 
R = { r k=0 , 

arcsinr k = 1 

generates an isometry along these straight lines and 
therefore there is a conserved quantity P of geodesic mo- 
tion attached to it, the specific linear momentum of the 
observer, 



±P 



d R = a 2 (t)f(r)f, 



(3) 



where the • denotes the inner product defined by the 
metric ([TJ) . The double sign is introduced in order to 
keep P positive. 

We just need another equation for i to complete the 
set and we may obtain it without resorting to Christoffcl 
symbols for the metric by using the unitarity condition 
©, 

5 = i 2 -a 2 (t)f 2 (r)r 2 . 



Restricting to future-pointing geodesies, t > (past- 
pointing geodesies are treated in a similar fashion), the 
whole set of geodesic equations is reduced to 



±- 



P 2 

p 



2 (t)f(r) 



(4a) 
(4b) 



Hence we see that there are basically three types of 
causal geodesies: radial lightlike ({ = 0, P / 0) and 
timelike geodesies (8 = 1, P ^ 0) and the comoving con- 
gruence of fluid worldlines (6 = 1, P = 0), which provide 
little information about the geometry of spacetime, since 
for them t = r regardless of the possible singularities in 
the universe. 



III. SINGULARITIES AT INFINITY 

Since singularities along causal geodesies at a finite to 
were considered in detail in [Til ] , we focus now on infinite 
values of coordinate time t. 

Singularities may appear also at t = ±oo if there are 
observers that reach these events in finite proper time. 
Unfortunately, it is not always possible to perform power 
expansions of the scale factor centered in t = ±oo, as it 
is done in [^, [ll| for finite t, since there are physically 
reasonable spacetimes with oscillatory scale factors, for 
instance, anti-de Sitter universes, for which the limit of 
a(t) is not defined when t tends to infinity. 

However, the question of when t = =boo is reached 
by geodesic observers in finite proper time can be easily 
solved. 

For lightlike radial geodesies we have 



t = 



P 



a(t')dt' = P(t~t ), 



la 



and therefore lightlike geodesies reach t = oo in finite 
proper time if and only if the integral 



a(t') dt' 



(5) 



is finite for sufficiently large t. That is, if a(t) is an 
integrable function at infinity. 

Comoving fluid wordlines with P — need not be con- 
sidered, since they reach t = oo in infinite proper time. 

Finally, we have timelike radial geodesies. In this case, 
proper time may be written again in terms of an integral 
of a(t) using (|4"a[) . 



dt' 



= T-T , 



(6) 
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and therefore these geodesies reach t — oo in finite proper 
time if and only if the improper integral 



df 



For non-tilted spatial geodesies in a hypersurface t 



t y/l + P 2 /a 2 {t>) 



is convergent for sufficiently large t. 

Obviously this can only happen if a(t) tends to zero at 
infinity, but it is not a sufficient condition. Since we may 
bound 



df 



i r°° 

= <pj t <f)df 



the integral for timelike geodesies is convergent if the one 
for lightlike geodesies is. 

Furthermore, since for large t and a(t) tending to zero 



1 



a{f) la 3 (f) 



P 



2 P 3 



is a telescopic series, the integral for timelike geodesies 
converges if and only if the one for lightlike geodesies 
does. Hence, all radial geodesies have the same regularity 
pattern. 

The analysis for t = — oo is entirely similar and so we 
have focused on the t — oo case. 

Since in most models the scale factor a(t) behaves 
asymptotically as a power of coordinate time, we start 
considering scale factors which behave close to infinity as 

2 

a(t)~c\t\ v , c> 0, w = 1. 

61] 

The equation for lightlike geodesies (j4aj) may be inte- 
grated close to infinity, 



(l + Tj)Pl 1/(1+n) 



t 



o P{t-t )/c 



»7 = -l, 



and provides valuable information, since t diverges when 
r tends to infinity for r\ > —1, whereas t diverges at finite 
proper time To if f] < — 1. 

The latter cases are quite interesting, since at To the 
geodesic reaches t — oo in finite proper time. Therefore, 
lightlike geodesies range from t = — oo (t = 0) to t = tq 
(t = oo) and are incomplete towards the future. 

This is not the interesting case, since it involves models 
starting at a Big Rip at t = 0. But If we consider t = 
— oo, lightlike geodesies range from t = To (t = — oo) to 
t = oo (t — 0) and are incomplete towards the past. This 
is the usual range in the suggested phantom models. 

As it has been said, the same behavior appears for 
timelike radial geodesies, with the difference that these 
actually end up at the Big Rip at t = in a finite proper 
time 

Not only causal geodesies, but also spatial geodesies 
show this feature. 



to, 



i = 0^t = t Q , P = a(t ), 



±- 



"o(*o)/(r) 



R= ±- 



•so 



a(t ) ' 



proper distance s is essentially the radial coordinate P, 
corrected by the expansion factor, as expected. 
But for tilted spatial geodesies, 



t = 



I p2 

aHt) 



- I => s ~ s = 



df 



t ^P 2 /a 2 {t') -1 



and for a(t) ~ c\t\ v for large t, this integral converges to 
a finite value if and only if 77 < — 1. Hence the length of 
these tilted spatial geodesies is also finite, even though 
the radial coordinate r diverges. 



IV. CURVATURE SINGULARITIES 

However, at t = ±00 all curvature scalar polynomials 
vanish, since they decrease as t~ 2 and this suggests a 
sort of Minkowskian limit. Therefore there is no scalar 
polynomial curvature singularity there. A pathological 
feature named imprisoned incompleteness, which appears 
in spacetimes like Taub-NUT, where geodesies are incom- 
plete without singular scalars of curvature, is not feasi- 
ble, since the spacetime has a cosmic time (a function 
with timelike gradient everywhere, the coordinate time 
f) and is therefore causally stable [3, . We might sus- 
pect that geodesic incompleteness could simply point out 
that the spacetime is not fully covered with the coordi- 
nate patch (fTJ) and that therefore it could be extendible 
beyond t = ±00. 

This is the case, for instance, of Milne universe, corre- 
sponding to k = — 1, ait) = t in ([T]). A suitable coordi- 
nate transformation 



T = t^l + r 2 , R = rt, 

shows that this model is just the portion of Minkowski 
spacetime inside the null cone T = R and therefore the 
apparent singularity at t = is due just to the choice of 
coordinates. 

A similar feature exhibits de Sitter spacetime in the 
parametrization that is usually used for inflation, k = 0, 

a(t) = eV A / 3 * ; which fulfils condition ([5]) and there- 
fore its radial geodesies reach t = —00 in finite proper 
time. However, again in this case it is possible to ex- 
tend the spacetime to a larger one, k = 1, a(T) = 

a/3/ A cosh ( y/ A/3 T^j with another change of coordi- 
nates [lH and hence the singularity at t = —oo is only 
apparent. 

Or for Schwarzschild spacetime, which in 
Schwarzschild coordinates appears to be singular at the 
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horizon at r = 2M, whereas this coordinate singularity 
disappears on extending it with Eddington-Finkelstein 
(2(il ] or Kruskal 21} coordinates. 

However, the null value of the scale factor in that limit 
suggests a point as a limit in this case. 

If we compute the Ricci tensor component along the 
velocity of the geodesic, an exotic behavior appears. 

For a radial lightlike geodesic, 



t = 



r • 1 

u = r = ± 



P 



RijU l u 3 



2P l 



2P 2 i] 



2kP 2 



~ c 2 t 2( n +l) c 4 t 4r, ' 

we take a look at the first term of the zero component of 
the Ricci curvature, since it is present regardless of the 
value of fc, 



2P Z \— - 



2 V 



1 



2P 2 T] 

c 2 t 2 (v+i) - (^ + 1)2 (r-r ) 2 ' 



(8) 



and find out that the Ricci curvature diverges when t 
approaches ±00 (r tends to t ) for 77 < — 1. 

The result invoked for the singularity-free de Sitter 

spacetime, k = 0, a(t) = ev 7 ^ 3 * is also recovered since 
in this case the expression for the Ricci curvature along 
the geodesies is zero, 

a 12 + k = aa". 

The remaining solutions for this equation are other 
parametrizations of the de Sitter spacetime, k — 

1, a(t) = y/S/Acosh(jAj3tj-, k = -1, a(t) = 
^3/Asinh (^y^ A/3tj and the anti-de Sitter spacetime, 

k = — 1, a(t) = -\/3/A cos A/3tj , but for a choice of 

time coordinate origin. None of them are decreasing at 
t = ±00, so they do not affect our results. 

A similar analysis may be performed for radial timelike 
geodesies, 



"' = '= 1 ' ~J- 



r • 1 

u = r = ± 



fa 2 ' 



RijU l u 3 



3a" 
a 



2P l 



(9) 



since for scale factors with r\ < — 1 and large values of \t\ 
the dominant term is the P 2 -term, which is the same as 
for lightlike geodesies. 

Hence the Ricci curvature diverges along both families 
of radial geodesies on approaching t = ±00. 

Hence we are to conclude that universes with scale fac- 
tor a{t) ~ c\t\ v , rj < —1, for large values of \t\ have a p.p. 
curvature singularity (curvature singularity along a par- 
allelly transported basis) [l8[ at t = ±00, though the 
scalar polynomials of curvature are zero there. 



That is, there is an actual curvature singularity at 
t = ±00 for these models, which correspond to w € 
(—5/3,-1), which is reached by the observers in finite 
proper time. Considering just expanding models of this 
type, all radial observers would trace their geodesic paths 
from the initial singularity at t = — 00 to the Big Rip at 
t = in a finite lapse of proper time. 

This result does not contradict the Penrose diagrams 
for these models shown in [22j , since conformal diagrams 
provide no information about distances, just about an- 
gles, but introduces a difference between models with 
w € (—5/3, —1) and those with w < —5/3 as it is shown 
in Fig. 1. 



r=0 




FIG. 1: Conformal diagram for a model with w G (—5/3, — 1): 
Timelike radial geodesies like a have finite length, whereas 
lightlike geodesies like b axe infinite towards the future and 
timelike geodesies like c are infinite towards the past. 



Furthermore, we may check the strength of these cur- 
vature singularities, which might be relevant, since other 
types of singularities, such as sudden singularities [7| (sin- 
gularities II and IV in [T3|) were shown not to be strong 
enough to disrupt finite objects [HJ and even have been 
suggested to be consistent with observations [HI]. 

Definitions of singularities related to curvature and 
geodesies refer to ideal point observers. When finite ob- 
jects are considered, tidal forces are relevant and it is 
interesting to check if they may destroy the object. In 
this case the singularity is considered to be strong [l3|. 
This qualitative concept has been stated rigorously by 
several authors [3, [H, [H, [H . 

For instance, in Tipler's definition [l4j a curvature sin- 
gularity is strong if the volume spanned by three Jacobi 
fields refered to a orthormal basis parallely-transported 
along a causal geodesic tends to zero at the singularity. 
Krolak's definition [15j just requires that the derivative 
of this volume be negative. 

There are necessary and sufficient conditions for the 
appearance of strong singularities [26] , that become quite 
simple to implement in the case of FLRW spacetimes, 
since the Weyl tensor vanishes 

With Tipler's definition a lightlike geodesic of velocity 
u comes up a strong singularity at tq if and only if the 
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integral 



dr' 



dT"R ljU l u } 



(10) 



diverges as r tends to tq. 

And with Krolak's definition a lightlike geodesic of ve- 
locity u comes up a strong singularity at tq if and only 
if 



dr' ' RijU l u° 



(11) 



diverges as r tends to To. 

Since the Ricci curvature component ([7J diverges as 
l/(r — r ) 2 , the integral of this term provides a logarith- 
mic divergence at To for r/ < — 1 with Tipler's definition 
and an inverse power divergence with Krolak's definition. 
Therefore, lightlike geodesies meet a strong singularity at 
t = ±oo if and only if r\ < — 1. The contribution of the 
curvature term diverges even faster when present. 

The previous conditions on integrals of Ricci compo- 
nents become sufficient conditions on dealing with time- 
like geodesies. Since the behavior of Ricci curvature 
has been shown to be similar for both families of radial 
geodesies for large \t\ and 77 < — 1, we learn that also 
radial timelike geodesies meet a strong curvature singu- 
larity at t = ±00. 

We may easily extend this result to non-power law 
growth/decrease of the expansion factor: 

• For a(t) growing or decreasing as l/\t\ or slower, 
radial geodesies reach t = ±00 in infinite proper 
time. 

• For a(t) decreasing faster than l/\t\, radial 
geodesies reach t — ±00 in finite proper time and 
therefore there is an actual strong curvature singu- 
larity there, except for de Sitter spacetime. 

These two cases include all situations for which the 
scale factor has a well-defined limit t — » ±00. Oscillatory 
scale factors may be treated directly with condition ([5]). 



V. DISCUSSION 

So far we have shown that FLRW cosmological models 
for which a(t) decreases faster than l/\t\ for large values 
of \t\ show a strong curvature singularity for t — > ±00, ex- 
cept for de Sitter spacetime. This is the case of phantom 
models with w € (—5/3, —1), family that includes mod- 
els compatible with observations, since w is estimated to 
be slightly below minus one [17| . 

Since the implications of these results are related to 
the past of the models instead of their future, it might 
seem a pointless discussion, for phantom models are in- 
tended to describe the future of the universe from now 
on. In the past other fields such as dust, radiation and 
the cosmological constant would be dominant and would 



prevent the appearance of the exotic curvature singular- 
ities described here. 

However, even though phantom models are not rele- 
vant to study the past of the universe, there are still 
consequences that are applicable to our present universe. 

We may consider, for instance, the total duration of a 
universe filled with a phantom field as experienced by a 
free- falling observer ©, 



T = 



dt 



y/i + P 2 /a 2 (t) 

1/7) 



dt 



c) Jo \/T 

p \ l/v B 

c 



x n dx 



,,2,, 



j_ 1 

2ri> 2 



2// 



v < -1, 



(12) 



-(c/Py/'it, using the hy- 



by the change of variable x 
pergeometric function Beta. 

We already know that this expression for the time span 
is finite for -q < — 1, but it can be made as small as desired 
by taking arbitrary large values of the linear momentum 
of the observer P. There is no lower bound, nor upper 
bound, which we know it is infinite for non-radial ob- 
servers. 

Though the calculation has been made for the dura- 
tion of the universe from the initial singularity to the 
Big Rip, it is clear that this result is also valid for the 
time span from the coincidence moment when phantom 
fields become the dominant component of the content of 
the universe to the Big Rip. That is, non-accelerated 
observers may shorten the time span to the end of the 
universe at will by increasing their linear momentum P. 

This feature is exclusive of dark energy models with 
w e (—5/3,-1), since a negative exponent r\ is required 
in (fT2)) for the decreasing behavior of T. For models with 
w < —5/3, 77 € [—1,0) the integral (fT2")) is divergent, 
since t = — 00 is actually at infinity, but the decreasing 
behavior is also exhibited for finite intervals of time up 
to the Big Rip singularity, though the interest on these 
models is so far quite limited. 

Another issue is the character of the singularity. Since 
a(t) tends to zero at t = ±00 for models with w G 
(—5/3,-1), this might suggest a sort of Big Bang sin- 
gularity, though endowed with exotic features. How- 
ever, the sign of the Ricci curvature measured by causal 
geodesies ([71 prevents this interpretation, since in (|S]) 
we see that it is negative (non-focusing) for negative 77 in 
flat models. 

In fact, for radial geodesies, 



dR 
~dt 



R 



±- 



P 



±- 



1 



t a{t)y/P*+a*(t) a(t)' 



R c± i?o i 



1 t 1 -^ 



c 1 



'I 
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the radial coordinate diverges for large \t\ in models with 
i] < 1, and shows that geodesies are indeed not focusing 
and therefore geodesies diverge instead of converge. 

With all these results in mind, we may refine the usual 
classification of singular events in models according to the 
value of w by introducing this new information, bearing 
in mind that none of these models is valid for the whole 
life of the universe, just for a fraction of it: 

1. Events with r\ > 1, —1 < w < —1/3: Quintessence 
models with a Big Bang singularity at t = 0. 

2. Events with r/=l, w = -1/3, k = -I: Milne-like 
models which may have weak or strong singularities 
at t = [lOU. 

3. Events with < r\ < 1, w > —1/3: Classical mod- 
els (dust, radiation. . . ) with a Big Bang singularity 
at t = 0. 

4. Events with 77 = 0: The menagerie of models which 
are either regular (de Sitter, anti-de Sitter. . . ) or 
possess sudden, freeze, pressure, higher derivative 
singularities as described in 0, [l(| [Tl[ , which may 
be weak or strong. 



5. Events with — 1 < 77 < 0, w < —5/3: Phantom 
models ranging from t = — 00 to i = where they 
meet a Big Rip singularity. Scalar perturbations 
of these models have been shown to lead to high 
inhomogeneity, which may prevent the formation 
of the singularity (27[. 

6. Events with r\ < —1, w 6 (—5/3,-1): Phantom 
models with a p.p. curvature singularity at t = —00 
which affects just radial geodesies and a Big Rip 
singularity at t — which does not show up for 
lightlike geodesies. 
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